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' The effect of strong magnetic fields on the equation of state (EoS) for compact stars 

I described with density dependent relativistic hadronic models is studied. A comparison with 

other mean-field relativistic models is done. It is shown that the largest differences between 
models occur for low densities and that the magnetic field affects the crust properties of star, 
namely its extension. 
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Cf^ ■ The study of very asymmetric nuclear matter is presently an important issue due to the ra- 



OO 
O 



INTRODUCTION 



dioactive beams which will be operating in the near future and which will allow the investigation 
of a region of the nuclear matter phase space unaccessible till recently. Asymmetric nuclear matter 
is of particular interest for the description of stellar matter of compact stars. 
^ ' Compact star properties depend a lot on the model used to describe the hadronic equation of 

state (EoS). In particular relativistic nuclear mean-field models [1, 2] are very popular to describe 
stellar matter because causality will always be satisfied. The imposition of constraints, both coming 
from measured star properties or from relativistic heavy ion collisions in the laboratory, is essential 
to test the different models [3]. 

Magnetars are neutron stars which may have surface magnetic fields larger that 10^^ G [4-6] and 
which were discovered at the x-ray and 7-ray energies (for a review refer [7]). They are identified 
with the anomalous x-ray pulsars (AXP) and soft 7-ray repeaters. Taking as reference the critical 
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field at which the electron cyclotron energy is equal to the electron mass = 4.414 x 10^ G 
we define B* = B/B^. It has been shown by several authors that the magnetic fields larger than 
B* = 10^ will affect the EoS of compact stars [8, 9]. In particular field-theoretical descriptions based 
on the non-linear Walecka model (NLWM) [1] were used and several parametrisations compared 
and it was shown that they had an overall similar behaviour. Very strong magnetic fields can only 
occur in very young compact stars before the magnetic field has decayed. Recently it was shown 
[10] from a 2D calculation of the cooling of magnetized stars that the magnetics fields and Joule 
heating have an important effect of maintaining compact stars warm for a longer time. This kind of 
simulations need the EoS of the crust. It is, therefore, important to make a study that shows when 
should the magnetic field be taken into account explicitly in the EoS of the crust. An unstable 
region in a wider density range will correspond to a larger crust and the properties of the star 
depending on the crust will be affected. It should, however, be pointed out the estimated surface 
magnetic fields of detected magnetars is obtained assuming that the lost of angular momentum is 
entirely due to dipolar radiation of magnetars, and until recently, the strongest estimated magnetic 
field is of the order of B* = 10^ and was detected in a quite young star, SGR 1806-20 [11]. 

At low densities relativistic models with constant coupling parameters have different behaviour 
from density dependent relativistic hadron models (DDRH) [12, 13]. These models have density 
dependent coupling parameters and have originally been parametrised so that relativistic Dirac 
Brueckner Hartree-Fock (DBHF) calculations for nuclear matter were reproduced [14]. Within 
DDRH models the symmetry energy does not increase linearly with density as parametrizations of 
the NLWM, like the ones studied in [8], and show a behavior closer to non-relativistic models, either 
some of the recent parametrizations of the Skyrme interaction like SLy230a [15] or NRAPR [16] or 
variational microscopic calculations [17]. Since nuclear matter is composed of two different fluids, 
namely protons and neutrons, the liquid gas phase transition can lead to an isospin distillation 
phenomenon which has been confirmed experimentally [18]. In [19, 20] it was shown that the 
distillation effect was described both by NLWM and DDRH models, but DDRH models did not 
predict an effect so strong as the first ones. Moreover, it was shown that DDRH models have at 
subsaturation densities a behaviour similar to non-relativistic models, namely models with Skyrme 
forces [21, 22]. We would like to test these models under different conditions namely /3-cquilibrium 
matter under strong magnetic fields. Stellar matter, as found in compact stars, under strong 
magnetic fields has already been studied before by many authors [8, 9, 23, 24]. 

The authors of Ref. [25, 26] have stressed the importance of including the scalar isovector vir- 
tual (5(ao(980)) field in hadronic effective field theories when asymmetric nuclear matter is studied. 
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Its presence introduces in the isovector channel the structure of relativistic interactions, where 
a balance between a scalar (attractive) and a vector (repulsive) potential exists. The introduc- 
tion of the S meson mainly affects the behaviour of the system at high densities, when, due to 
Lorentz contraction, its contribution is reduced, leading to a harder EoS at densities larger than 
~ 1.5 po [26]. In [27] the effect of this meson on the properties of compact stars were studied and it 
was shown that the EoS of hadronic matter would become stiffer with its presence. In relativistic 
models there is a proton-neutron mass splitting only if the scalar isovector (5-meson is included. 
This occurs for the DDRH(5 parametrization we consider. For this model M* < M* in neutron rich 
nuclear matter [28] . A similar behavior is predicted by the Skyrme interaction SLy230a but an op- 
posite behavior is obtained with other parametrizations of the Skyrme interaction [29] . In a recent 
work [30] the behavior of the proton-neutron mass splitting in different relativistic nuclear models 
was analysed and it was shown that some point-coupling models without mesons [31] predict larger 
neutron masses in neutron rich matter. The proton- neutron mass splitting is a present topic of 
discussion and the forecoming experiments with radioactive beams will allow the clarification of 
this point. We will investigate whether this mass spitting has some effect on stellar matter with 
strong magnetic fields. 

In the present paper we will study the behaviour of stellar matter described within DDRH, both 
with [32] and without [13] the isovector-scalar (5-meson, under very strong magnetic fields. The 
results are compared with previously studied models, the parametrisations GM3 [33] and TMl [34] 
of NLWM. These two models including strong magnetic fields, have been discussed in [8] and [23]. 

In section II we make a brief review of the model and EoS under the effect of a magnetic field. 
In section III the formalism is generalised to include the (J-meson. Results are discussed in section 
IV and conclusions are drawn in section V. 



For the description of the EoS of neutron star matter, we employ a field-theoretical approach 
in which the baryons (neutrons, n, and protons, p) interact via the exchange oi a — uj — p mesons 
in the presence of a uniform magnetic field B along the z-axis. The Lagrangian density of the 
relativistic TW model [12, 13] can be written as 



The baryons {b=n, p), leptons (/=e, p), and mesons {a, lo and p) Lagrangians are given by 



II. THE FORMALISM 




(1) 



b=n,p 



l=e,ix 
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Ly, = (^ilf^d'' - qbJf^Af" -mb + T^a - T^Tm^'' - ^'^pTsblt^p'" - \l^Nnb(J^,uF^'''^ *b 
>Cz = i^i {il^id^ - qnt^A'' -mi)^i 

- \f>^^F^, + \mlp,(f - \P^^P^. (2) 

where ^6 and if^i are the baryon and lepton Dirac fields, respectively. The nucleon mass and isospin 
projection for the proton and neutrons arc denoted by rrib and T^b = il, respectively. The mesonic 
and electromagnetic field strength tensors are given by their usual expressions: 0^,^ = d^uj^ — dyUn, 
P^v = d^j^pi, — dupn, and F^,^ = d^Ai, — duA^. The nucleon anomalous magnetic moments are 
introduced via the coupling of the baryons to the electromagnetic field tensor with cr^i^ = | [7/i, 7i/] 
and strength Kb with k„ = —1.91315 for the neutron and Hp = 1.79285 for the proton, respectively. 
The electromagnetic field is assumed to be externally generated (and thus has no associated field 
equation), and only frozen-field configurations will be considered. The density dependent strong 
interaction couplings are denoted by F, the electromagnetic couplings by q and the nucleon, mesons 
and leptons masses by m. The parameters of the model are the nucleon mass nib = 939 MeV, 
the masses of mesons m^, m^;, rup and the density dependent coupling parameters which are 
adjusted in order to reproduce some of the nuclear matter bulk properties and relations with 
DBHF calculations [14] , using the following parametrisation 

^i{p) = 'i^i{Psat)fi{x), i = a,(jj (3) 



with 



where x = p/psat and 



l + bi{x + d, 



\2 



fi (^) = «i 1 ^ / ^^^2 (4) 
l + Ci{x + di) 



^p{p) = ^piPsat) exp [-ap{x - 1)] , (5) 

with the values of the parameters rrii, Fj, Oj, bi, Ci and di, i = a,u;, p given in [13]. Other possibilities 
for these parameters are also found in the literature [35]. 

The field equations of motion follow from the Euler-Lagrange equations. Prom the Lagrangian 
density in Eq. (1), we obtain the following meson field equations in the mean-field approximation 

ml {a) = T„{p; + pi)=T^p' (6) 
m2(a;°) = [p; + p^ = T^pb (7) 
,{p') = lrp{p;-pl) = lT,ps (8) 



ml 
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and the Dirac equations for nucleons and leptons are given by 



{i-f^d^" - qa^A^" - mi) il^i = 



(9) 
(10) 



where the effective baryon masses are given by 



mt = mb- T„a 



(11) 



and is the scalar number density. In charge-neutral, /3-equilibrated matter, the conditions 



(12) 



and 



(13) 



should be satisfied. 

The energy spectra for protons, neutrons and leptons (electrons and muons) are given by 

El,s = + ( '^^QpB - sfXNKpB) ' + r^a;° + ^r^p^ + (14) 

= sjkl + ^kl + kl- siXNi^nB) ' + r^a;° - ^r^p^ + (15) 

Ei^, = ^Jkl + rrq + 2u\qi\B (16) 

where v = n + ^ — sign{q)^ = 0,1,2,... enumerates the Landau levels of the fermions with 
electric charge q, the quantum number s is +1 for spin up and —1 for spin down cases, and the 
rearrangement term is given by 



^° " dp 



dp 



p a. 



(17) 



The expressions of the scalar and vector densities for protons and neutrons are given by [8] 



qpBm*p yjrnf + 2uqpB - sfiNKpB 



v=0 s 



m*^ + 2uqpB 



Pn 



47r2 
qpB 



-mS In 



s ""n 



rrip + 2vqpB — spjyUpB 



1^=0 s 
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Pn 



27r2 



E 



[^F,s) ^ ~ i^siiNHnB ( mnkp^s + E'l? ( arcsin ( ^ 



and the vector densities for leptons are given by 



Pi 



\B 



27r2 



(18) 



(19) 



:/=0 s 

where k^p^g-, kp^ and /c^^^, arc the Fermi momenta of protons, neutrons and leptons, which are 
related to the Fermi energies E'^p, E^ and E^p as 

n2 



LP2 _ pp2 



kf, = Ef-ml 



with 



k%,s = e3 - {mf + 2iy\qi\B) , I = e, 



(20) 



(21) 



The summation in v in the above expressions terminates at Vmax, the largest value of u for which 
the square of Fermi momenta of the particle is still positive and which corresponds to the closest 
integer from below defined by the ratio 



2\qi\B 



leptons 



{E^p + s /ijv Bf - m*p 



*2 



2\qp\B 

The chemical potentials of baryons and leptons are defined as 



, protons. 



fib = E''p + r^a;° + ^Tpmp"^ + E, 



fil = E^p = Jkf^^ + mf + 2v\qi\B. 



(22) 
(23) 



We solve the coupled Eqs. (6)-(13) self-consistently at a given baryon density in the presence of 
strong magnetic fields. The energy density of neutron star matter is given by (the index "m" refers 
to matter) 

^ ^ 111 

£m = £6 + X] + + 2"'^'^o + 2^y^ ^^^^ 

b=p,n l=e,ij, 

where the energy densities of nucleons and leptons have the following forms 



i>=0 s 



f^F,u,sEF + {\J^*p + 2z^gp-B - SIlNKpB^ In 



m*2 + 2vqpB — sf^Ni^pB 
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rrir 



^F,s + -^f 



I-B 



47r2 



EE 
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m? + 2i'\qi\B 



(25) 



The pressure of the system is obtained from the expression 



(26) 



where the charge neutrahty and /3-equihbrium conditions arc used to get the last equahty. Note 
that the contribution from electromagnetic fields to the energy density and pressure, £f = Pf = 
should be taken into account in the calculation of the EoS. 



III. INCLUDING ISOVECTOR-SCALAR MESONS 



To investigate the influence of the (5-meson we have included in the TW model the isovector- 
scalar meson term [32], which has density dependent coupling parameters. The Lagrangian density 
reads 



l=e,iJ, 



(27) 



where the baryon {b=n, p), lepton {l=e, /x), and meson {a, to, p and 6) Lagrangian are given by 

A = ^{d^6d>^6-mj?), (28) 

with Cb and Cm defined in Eq (2). Ts and ms are, respectively, the coupling constant of the 
6 meson with the nucleons and its mass. For Tu and we take the parametrisations given in 
Eqs. (3) and (4). For Tp and Ts, we use the parametrisation [36] 

P 



Ti{p) = Ti{psat)fi{x), X = 



Psat 



with 



fi{x) = Ui exp [-bi{x - 1)] - Ci (x - di) , i = p,S 



(29) 



(30) 



and the parameters Oj, bi, ci and di defined in Table I. 
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i Ti 


ai 


bi 




p 11.727 


0.095268 


2.171 


0.05336 17.8431 


5 7.58963 


0.01984 


3.4732 


-0.0908 -9.811 



TABLE I: Parameters of the DDRH5 model. 

From the Lagrangian density in Eq. (27), we obtain the meson field equations (6), (7), (8) plus 
an equation for the c^-meson in the mean-field approximation 

r4{5^) = Vs{pl-pt,)=Tspl (31) 

and the Dirac equations for nucleons are given by 

{ilixd^ - QblixA^ - {mi, - r^a - TsTsbSs) - r^7ow° 

-^rpr367o/ - 70^0 - l^^Nl^b<T^c.F''''^ = 0. (32) 

The effective baryon masses, in this case, are given by 

ml = mh-T^a-TsbTs53, (33) 

and differ for protons and neutrons. In charge-neutral, /3-equilibrated matter, the conditions 
Eq. (12) and Eq. (13) apply. 

The energy spectra for protons, neutrons and leptons are given by Eqs. (14), (15) and (16) with 
the rearrangement term, in this case, given by 

So = -Q^PbO^o + ^P3y - ^ + -Q^PsSs. (34) 

The expressions of the scalar and vector densities for protons and neutrons, Fermi momenta, 
chemical potentials still hold. The energy density of neutron star matter is, now, given by 

£m= ^ £b+^ £1 + + \mi<A + ^rn^pl + ^m^5| (35) 

b=p,n l=e,ix 

including an extra term for the (5-meson. For the pressure Eq. (26) holds. 

IV. RESULTS AND DISCUSSION 

In the present section we discuss the EoS of stellar matter obtained within TW and DDRH(5 and 
compare them with other previously studied models GM3 and TMl. We will pay a special attention 
to the behaviour of the EoS at subsaturation densities in order to understand how could strong 
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magnetic fields affect the crust of compact stars by extending or reducing tlie non-liomogeneous 
phase. 

The properties of the isovector channel of the models have an important role on the properties 
of the EoS of very asymmetric matter and therefore we plot, for reference, in Fig. 1 the symmetry 
energy of the four models under study. It is seen that the models have quite different behaviours at 
large densities and this will reflect itself on the EoS. In fact the largest proton fractions occur for 
the models with the largest symmetry energy since a very asymmetric system will be energetically 
defavoured. While the DDRH(^ is the one with the lower symmetry energy at lower densities, 
around two times the saturation density its symmetry energy crosses the corresponding curve for 
the TW. This is due to the saturation of the 5 field as discussed in [26]. We may expect that this 
effect will have influence on the properties of the EoS. GM3 and TMl have quite high symmetry 
energy esym which originates a large proton fraction at subsaturation densities and allows for the 
direct URCA process. We will see, however that the proton fraction is determined also by the 




FIG. 1: (Color online) The symmetry energy of the models under study. 
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FIG. 2: (Color online) EoS for stellar matter without (left) and with (right) the nucleon anomalous magnetic 
moment. Several models are considered. 

In Fig. 2 we compare the EoS of TMl, GM3, TW and DDRHJ for several magnetic field 
intensities. As discussed in [8] the magnetic field makes the EoS softer when the anomalous 
magnetic moment (AMM) of the nucleons is not considered. For B* = 10^ the strongest effect 
occurs at low densities typical of the crust of the star. It is therefore important to investigate the 
influence of the magnetic field on the crust properties. For S = it is shown that DDRHJ presents 
a small liquid-gas phase transition. This is certainly due to the small symmetry energy this model 
has at low densities. From a qualitative point of view all models behave similarly, also when the 
AMM is included. In this case the EoS becomes stiffer as the magnetic field increases except for 
the small densities. Inparticular the curve for B* = 10^ gets softer than the B = EoS for energy 
densities smaller than 70 MeV/fm^, when the effect of the Landau quantisation is stronger than 
the AMM contribution. 

In order to have a more quantitative comparison, in Fig. 3 we plot for each magnetic field the 
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FIG. 3: (Color online) EoS for stellar matter without (left) and with (right) the nucleon anomalous magnetic 
moment. Several magnetic fields are considered. 

four models. If the anomaly is not considered the largest differences occur for the lower densities: 
at i? = TMl and GM3 are the softest at low densities and become the hardest at high densities; 
DDRH(5 is the softest EoS at intermediate energies. The B = behaviour is determined by the 
symmetry energy, and the trend of the EoS follows the relative behaviour of the symmetry energy 
curves. However the relative stiffness of the EoS depends on the intensity of B. Contrary to the 
lower values of B*,0 and 10^, for the highest value GM3 is stiffer than TMl and for B* = 5 x 10^ 
both EoS have very similar behaviours. Although DDRH5 EoS is softer than the TW EoS below 
e = 200 MeV/fm^ for small and large magnetic fields, for B* = 10^ and 5 x 10^ it becomes harder 
or similar. If the anomaly is included the differences between the models are much smaller mainly 
for B* < 10^. For these values, models coincide at intermediate and high densities. For the low 
densities, differences arise for B* ~ 10^ or lar ger. 

In Fig. 4 only the region p < l.bpo is plotted. The magnetic field increases the binding at 
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FIG. 4: (Color online) EoS at low densities for stellar matter without (left) and with (right) the nucleon 
anomalous magnetic moment. Several magnetic fields are considered. 

low densities and a negative pressure may still occur beyond 2po- The most bound matter occurs 
for DDRH5 and this is due to the low symmetry energy this model has at intermediate densities. 
This means that the crust extends itself to higher densities and a higher value of the density 
will characterise the inner edge of the crust. The inclusion of the AMM reduces this effect but 
binding still occurs and there are also regions of negative compressibility and/or pressure defining 
an unstable region. It is clearly seen that the models behave differently in these range of densities 
when the B field increases. GM3 becomes harder than TMl for B* > 5 x 10^. Also the relative 
behavior of TW and DDRH(5 change. 

In order to better understand the low density behavior we have determined the upper density 
limit of the mechanical instability, the density at which the incompressibility becomes zero, as a 
function of the magnetic field. Fig. 5a), and calculated the associated proton fraction Fig. 5b). For 
/3-equilibrium nuclear matter with no magnetic field there are no mechanical instabilities even if 
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FIG. 5: (Color online) Low density mechanical instability. Bottom: The upper limit of the mechanical 
instability for different magnetic fields as a function of density; top: the corresponding proton fraction; left: 
no AMM; right: including AMM. 

we only consider neutron-proton (np) matter and do not add the electron contribution, for which 
the incompressibility is always positive. In Table II we give the densities for which the EoS for 
/3-equilibrium matter at a zero magnetic field crosses the thermodynamical spinodal for np matter. 
As shown in [37] these numbers give an order of magnitude of the upper limit of the transition 
density: neither the Coulomb force nor the finite range of the nuclear force is taken into account. 
They are slightly larger than the corresponding values obtained from the crossing of the dynamical 
spinodal with the EOS for neutron-proton-electron (npe) matter, which define a lower limit, and 
of the same order of magnitude of the results obtained from the transition of a pasta phase to a 
homogeneous phase [37]. 

From Fig. 5 we see that the magnetic field, if strong enough, may change this picture, more 
strongly if the AMM is not considered. In this figure we show both the mechanical instability upper 
bound for np matter and also for npe matter. Of course in the last case the instability region is 
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TABLE II: Predicted density and pressure at the inner edge of the crust of a compact star at zero tem- 
perature, as defined by the crossing between the thermodynamical instability region of np matter and the 
/3-equiUbrium condition for homogeneous, neutrino-free stellar matter. 





Pfc(fm-3) 


P™(McVfm"^) 


TMl 


0.069509 


0.50288 


GM3 


0.068762 


0.35644 


TW 


0.084955 


0.52246 


DDRH(5 0.085038 


0.12855 



smaller due to the high incompressibility of the electron gas and occurs at larger proton fraction 
when the symmetry term contribution is smaller. The existence of a mechanical instability region 
in the presence of the magnetic field has two reasons: a) due to the existence of Landau levels 
the nucleonic pressure does not increase so fast with density and b) at large magnetic fields the 
proton fraction increases and the symmetry repulsive term in the energy density is not so strong. 
The model DDRH(5 is the one showing the largest instability ranges and at least, for the no AMM 
calculation, it also predicts the smaller proton fractions. However, when the AMM contribution 
is introduced the proton fraction of all models behave in a very similar way except for the larger 
fields. In fact we should perform a calculation which includes the Coulomb interaction and surface 
energy, but according to studies done in [21, 38, 39], the spinodal which includes these effects would 
be larger than the mechanical instability region we have calculated. A complete study of the low 
density region, namely the spinodal surface that limits the non-homogeneous phase, needs to be 
done. 

The fraction of protons and muons for the models and field intensities discussed are given in 
Fig. 6 as a function of density. For comparison we also include the B = results. The effect of the 
magnetic field is not very large for B* = 10^, except at low densities, p < 2pQ, when protons are 
totally polarised. In Fig. 7a) we show the occupied Landau levels as a function of density and in 
Fig. 7b) the neutron polarization for the calculation including AMM. For B* = 10^ {B* = 5 x 10^) 
the second Landau level starts being occupied only for densities above ~ 2pQ (~ 8po). GM3 has a 
similar behaviour to TMl and DDRH(5 an intermediate behaviour between TMl and TW. From 
Fig. 7b) we conclude that while at B* = 10^ neutrons are only slightly polarized at B* = 10^ they 
are totally polarized for densities below 8po . This total neutron polarization favours an increase of 
the proton fraction when AMM is included. 

The fraction of protons within the different models is determined by the symmetry energy of 
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FIG. 6: (Color online) Proton and muon fraction for stellar matter without (left) and with (right) the 
nucleon anomalous magnetic moment. Several magnetic fields are considered. 

the respective model for the density considered. This explains why for B* < 10^, DDRHd has the 
smallest fraction for p < 2.5 pQ. For p > 2.5po it is TW which has the smallest symmetry energy 
and the smallest proton fraction. For high fields, B* > 10^ the fraction of protons is larger than the 
fraction of neutrons and again the symmetry energy defines the models with the largest fraction: 
DDRH5 for p < 3po and TW for densities larger than 3po- At B* = 5 x 10^ all proton fraction lie 
between 0.4 and 0.5 for densities below 4/Jo- The relative fraction of protons is then determined by 
the effective mass: GM3 has the lower fraction due to its larger mass. On the other hand DDRH5 
has the smallest proton mass a largest proton fraction. This effect is even larger for B* = 10^. 
The inclusion of the AMM has important effects for B* > 10^. For densities p < 2pQ the behaviour 
of the different models is still distinguishable but for larger densities the behaviour of all models 
is mostly determined by the magnetic field intensity. The magnetic field gives rise to an onset of 
muons at larger densities when the anomaly is not included because the magnetic fields favours a 
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FIG. 7: (Color online) Density dependence of the proton and neutron polarization for TMl and TW and 
different values of the magnetic field: a) Landau levels occupied by protons without (left) with (right) AMM. 
The dashed lines indicate the onset of spin down for each Landau level in the calculation with AMM; b) 
neutron polarization, full (dashed) lines correspond to spin down (up). 

larger proton fraction and therefore a smaller electron chemical potential. With the introduction 
of the AMM, the electron chemical potential increases slightly which explains the onset of muons 
at lower densities. 

In Fig. 8 the effective mass of the nucleons within the models under study are shown as a 
function of the density for different intensities of the magnetic field, with and without the AMM. 
As discussed in [8] when the AMM is not taken into account and the magnetic field becomes more 
intense, the effective mass reduces faster with an increase of the density. For B* = 10^ the effect 
of the magnetic field is still small, but for larger values it becomes more important. For B* = 10^ 
there is total polarisation only for p < 2po while for B* = 10® there is total polarisation for the 
all range of densities considered. The GM3 model has the largest mass for all fields considered. 
At saturation and without magnetic field this model predicts an effective mass of 0.78 M while 
for all the others the effective mass at saturation is 0.6 M or smaller. The large values of the 
effective mass within GM3 justify the smaller proton fractions for the most intense magnetic fields. 
On the order hand DDRH(5 shows the fastest reduction of the effective mass with density. For 
DDRH5 we show both the proton and the neutron mass. In this model, the mass of the most 
abundant nucleon, neutron for B* = 10^, 0.5 x 10® and proton for B* = 10®, or larger, behaves like 
the nucleon mass in all the other models (except GM3), having slightly smaller values. The mass 
of the less abundant nucleon is quite higher, similar to the nucleon mass in GM3 for the smaller 
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FIG. 8: (Color online) Nucleon effective mass for stellar matter without (left) and with (right) the nucleon 
anomalous magnetic moment. 



magnetic fields. If the AMM is included the effective masses saturate quite fast at a non-zero value 
which corresponds to 771^ — sfiNf^iB, i = p,n. Models with density dependent couplings saturate 
faster than the other two. 

The behaviour of the nucleon masses is better understood from the behaviour of the scalar field 
with density. In Fig. 9 we plot the a and 6 fields, more precisely Tua and — F^^a as a function of 
the density for several magnetic field intensities. If the AMM is not included the (5-field changes 
sign for the two most intensive fields considered. This reflects the existence of a larger fraction of 
protons than neutrons. The inclusion of AMM reduces this effect and only for a restricted range 
of densities below 3 — ApQ. This is due to the larger effective mass of the protons. We also see 
that the a field increases faster with density for the larger magnetic fields, giving rise to a faster 
saturation of the effective mass. If the AMM is taken into account the saturation of a field does 
not occur at the baryonic mass value but at the value of the renormalised baryonic mass by the 
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FIG. 9: (Color online) The cr and 5-fields for stellar matter without (left) and with (right) the nucleon 
anomalous magnetic moment. 

AMM, see Eq. (21). 

V. CONCLUSIONS AND OUTLOOKS 

In the present paper we have compared the EoS for stellar matter made out of protons, neu- 
trons, electrons and muons in the presence of very strong magnetic fields. In particular, we have 
studied the EoS obtained within DDRH models with and without the scalar-isovector meson 6 and 
compared with other models previously studied: GM3 and TMl. 

It has been shown that, although the overall behaviour of all the models is similar, at low 
densities, p < 3po, the models show the largest differences. In particular, it has been shown that 
the low density instability region increases a lot as the magnetic field increases when the AMM 
is not taken into account. Although the AMM reduces this behaviour there remains a region of 
instability not present for the magnetic field free matter. The larger range of instability is partially 
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due to the larger proton fractions and partially due to the appearance of Landau levels. 

Moreover DDRH5 model shows a low density instability, which will give rise to a liquid-gas 
phase transition, even for B = 0. This could be due to the low symmetry energy it has at 
low and intermediate densities. It is also the DDRH5 model that shows the largest changes of 
behaviour with the magnetic field and the density, both due to the presence of the 5-meson and 
the density dependence of the coupling parameters. This is particularly clear with the EoS and the 
proton and muon fractions. In fact, the DDRH models are the ones that predict smaller (higher) 
proton fractions at intermediate densities for B* < 10^ {B* > 10^), reflecting the behaviour of 
the symmetry energy, the proton effective mass with density and the difference between proton 
and neutron masses. Although the EoS of the different models do not differ so much, properties 
of the star sensitive to the proton fraction will distinguish the different models. These may be 
the neutrino interaction with hadronic matter, it is larger for larger neutron fractions, or the 
pairing properties of stellar matter which affect neutrino emissivities and specific heat. The stars 
more sensitive to the differences between the models are the low mass ones with M~ I.OMq. In 
particular, we have shown that the proton-neutron mass spitting present in DDR(5 parametrization 
has noticeable effects on the proton fractions predicted by this model, changing at intermediate 
densities from the smallest ones to the largest ones. 

In the present study we have not considered strangeness. At large densities about two times 
the saturation density we may expect the onset of hyperons [2], or kaon condensation [40]. The 
effect of strong magnetic fields on the onset of strangeness has been discussed in [41, 42] and it 
was shown that both the hyperon or the kaon condensate onset occurs at larger densities, > 5po, 
in the presence of strong fields. So we may consider in the present discussion that results will not 
be affected by the strangeness degree of freedom below 5 po- It remains, however to be checked if 
the density dependence of the baryon-meson couplings will have an effect on the hyperon onset. 

The study of the effect of the magnetic field on the low density instabilities is of particular 
interest: the way the clusterization occurs and the extension of the crust will affect the cooling 
and conduction properties of the star. It was shown in the present work that the properties of the 
crust under strong magnetic fields are sensitive to the EoS used. A detailed study of the effect of 
strong magnetic fields on the low density instabilities of nuclear matter is being carried on. 
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